Table I. Finite-iV values of e^ 1/2 , + )(N)/J, e ( i /2 _)(iV)/J, e^ 1/2 ,+y(N)/J, and e^/ 2 ,+)(N)/J for d = 0.0, 
of J'/J. Note that e (1/2 , + y(N) =e (3/2 ,+)(N) when J'>0 and e {3/2j+) (N) = e (1 / 2)+ )(iV) when J'<0. 



(a) e (1/2i+) (iV)/J 



J'/J 

I 


iV = 5 


N = 7 


iV = 9 


iV = ll 


iV = 13 


iV = 


-0.30 


-6.054182 


-8.848086 


-11.640230 


-14.433542 


-17.228587 


-20.0 


-0.20 


-5.974922 


-8.771783 


-11.565739 


-14.360276 


-17.156197 


-19.9 


-0.10 


-5.899732 


-8.699949 


-11.495923 


-14.291798 


-17.088660 


-19.8 


0.00 


-5.830213 


-8.634532 


-11.432932 


-14.230359 


-17.028266 


-19.8 


0.10 


-5.963368 


-8.759273 


-11.552826 


-14.347250 


-17.143216 


-19.9 


0.20 


-6.098600 


-8.887590 


-11.677328 


-14.469455 


-17.263959 


-20.0 


0.30 


-6.235457 


-9.018623 


-11.805252 


-14.595543 


-17.388888 


-20.1 



(b) e (1/2 ,_)(7V)/J 



J'/J 


iV = 5 


N = 7 


iV = 9 


iV = ll 


iV = 13 


iV = 


-0.30 


-5.309153 


-8.328417 


-11.246685 


-14.117088 


-16.962018 


-19.7 


-0.20 


-5.295485 


-8.315081 


-11.232456 


-14.101753 


-16.945633 


-19.7 


-0.10 


-5.286686 


-8.306572 


-11.223408 


-14.092010 


-16.935219 


-19.7 


0.00 


-5.283567 


-8.303576 


-11.220229 


-14.088587 


-16.931557 


-19.7 


0.10 


-5.286954 


-8.306811 


-11.223658 


-14.092287 


-16.935527 


-19.7 


0.20 


-5.297593 


-8.316962 


-11.234436 


-14.103943 


-16.948071 


-19.7 


0.30 


-5.316038 


-8.334613 


-11.253235 


-14.124347 


-16.970116 


-19.8 



(c) 6 {1/2 , + y(N)/J 



J'/J N = 5 N = 7 N = 9 N = ll N = 13 N = 



-0.30 -5.451354 -8.298314 -11.125222 -13.942842 -16.755529 -19.5 

-0.20 -5.573010 -8.401564 -11.215359 -14.023132 -16.828188 -19.6 

-0.10 -5.699777 -8.514645 -11.319475 -14.121061 -16.921741 -19.7 

0.00 -5.830213 -8.634532 -11.432932 -14.230359 -17.028266 -19.8 



(d) 6 (3/2 , +) (iV)/J 

J'/J N = b N = 7 N = 9 N = ll N=13 N = 



0.00 -5.830213 -8.634532 -11.432932 -14.230359 -17.028266 -19.85 

0.10 -5.768795 -8.578704 -11.380467 -14.180017 -16.979284 -19.7/ 

0.20 -5.719142 -8.537588 -11.344937 -14.148268 -16.950089 -19.7E 

0.30 -5.686463 -8.518738 -11.336204 -14.147241 -16.955021 — 19.7C 



Table II. Estimated values of A/^ A (i/2,+)'' and A( 3 / 2 ,+) for d = 0.0 and d = 
Note that A( 3 / 2 ,+) = ^(1/2,+y when d = 0.0. 



d=0.0 



d = 0.2 



J'/J 



A (l/2,-) A (l/2,+)' A (l/2,-) A (l 



(1/2,+)' 



k (3/2,+) 



-0.20 


-0.2015 


-0.3045 


-0.1728 


-0.2929 


-0.2609 


-0.15 


-0.1560 


-0.2337 


-0.1379 


-0.2307 


-0.2056 


-0.10 


-0.1061 


-0.1580 


-0.0944 


-0.1585 


-0.1409 


-0.08 


-0.0855 


-0.1269 


-0.0756 


-0.1277 


-0.1134 


-0.06 


-0.0646 


-0.0954 


-0.0561 


-0.0957 


-0.0852 


-0.04 


-0.0435 


-0.0637 


-0.0361 


-0.0623 




-0.02 


-0.0223 


-0.0319 








0.02 


0.0202 


0.0319 








0.03 






0.0309 






0.04 


0.0415 


0.0636 


0.0434 


0.0607 


0.0549 


0.06 


0.0626 


0.0952 


0.0659 


0.0953 


0.0836 


0.08 


0.0836 


0.1265 


0.0869 


0.1290 


0.1106 


0.10 


0.1044 


0.1573 


0.1072 


0.1623 


0.1365 


0.15 


0.1555 


0.2314 


0.1565 


0.2434 


0.1959 


0.20 


0.2048 


0.2990 


0.2036 


0.3179 


0.2444 



Table III. Estimated values of the coefficients cl(m,p) an d ^(m,p)- 
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bjQ' 
< 



d 


a (l/2,-) 


6 (l/2,-) 


0(1/2,+)' 


6(1/2,+)' 


a (3/2,+) 


& (3/2,+) 


0.0 


1.065 


0.018 


1.599 


-0.063 


1.599 


-0.063 


0.1 


1.036 


0.217 


1.604 


0.163 


1.539 


-0.290 


0.2 


1.026 


0.418 


1.618 


0.293 


1.426 


-0.210 


0.3 


1.014 


0.470 


1.628 


0.361 


1.325 


-0.117 


0.4 


0.990 


0.436 


1.627 


0.467 


1.230 


-0.022 
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Table IV. Finite-iV values of e ( i /2)+) (iV)/J, e ( i /2 _)(iV)/J, e W2 ,+y(N)/J, and e^/ 2 ,+)(N)/J for d = 0.2, 
of J'/ J. 



(a) e (1/2)+) (iV)/J 



J'/ J 


N = 5 


N = 7 


N = 9 


N=ll 


N=13 


N = 


-0.30 


-5.475582 


-8.012057 


—10.546742 


-io Anno /~\ A 

—13.082394 


-15.619570 


-18.1 


-0.20 


-5.394812 


-7.933029 


-10.468686 


-13.004973 


-15.542612 


-18.0 


-0.10 


-5.320835 


-7.860148 


-10.396151 


-12.932553 


-15.470260 


-18.0 


0.00 


-5.275712 


-7.814208 


-10.348532 


-12.882930 


-15.418570 


-17.9 


0.10 


-5.346869 


-7.886106 


-10.423171 


-12.961113 


-15.500510 


-18.0 


0.20 


-5.471736 


-8.006562 


-10.541781 


-13.078873 


-15.617811 


-18.1 
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(b) £(i/2, 


-)l iV >/ J 












J'/ J 


N = 5 


N = 7 


N = 9 


N = ll 


iV=13 


N = 


-0.30 


-4.652876 


-7.417147 


-10.081405 


-12.698384 


-15.290182 


-17.8 


-0.20 


-4.639113 


-7.403702 


-10.067030 


-12.682847 


-15.273516 


-17.8 


-0.10 


-4.630245 


-7.395117 


-10.057884 


-12.672970 


-15.262916 


-17.8 


0.00 


-4.627100 


-7.392093 


-10.054669 


-12.669498 


-15.259185 


-17.8 


0.10 


-4.630519 


-7.395361 


-10.058140 


-12.673254 


-15.263234 


-17.8 


0.20 


-4.641266 


-7.405623 


-10.069054 


-12.685093 


-15.276033 


-17.8 


0.30 


-4.659906 


-7.423472 


-10.088095 


-12.705821 


-15.298533 


-17.8 



(c) e {1/2 , +y (N)/J 



J'/J 


N = 5 


N = 7 


N = 9 


N = ll 


iV = 13 


N = 


-0.30 


-4.813371 


-7.412625 


-9.991555 


-12.560630 


-15.124107 


-17.6 


-0.20 


-4.931077 


-7.511422 


-10.076178 


-12.633803 


-15.187583 


-17.7 


-0.10 


-5.050468 


-7.617882 


-10.174052 


-12.725463 


-15.274475 


-17.8 


0.00 


-5.149004 


-7.709056 


-10.261722 


-12.811685 


-15.360390 


-17.9 


0.10 


-5.139953 


-7.693920 


-10.240734 


-12.785313 


-15.329197 


-17.8 


0.20 


-5.087907 


-7.643786 


-10.191231 


-12.736204 


-15.280593 


-17.8 


0.30 


-5.042013 


-7.604356 


-10.155913 


-12.704125 


-15.251613 


-17.7 



(d) e (3A+ )(iV)/J 



J'/J 


N = b 


N = 7 


N = 9 


iV = ll 


iV = 13 


N = 


-0.30 


-5.355288 


-7.901662 


-10.446271 


-12.991136 


-15.536428 


-18.0 


-0.20 


-5.280562 


-7.830554 


-10.377293 


-12.923541 


-15.469803 


-18.0 


-0.10 


-5.211118 


-7.765469 


-10.314753 


-12.862610 


-15.409957 


-17.9 


0.00 


-5.149004 


-7.709056 


-10.261722 


-12.811685 


-15.360390 


-17.9 


0.10 


-5.097223 


-7.665435 


-10.223245 


-12.776551 


-15.327447 


-17.8 


0.20 


-5.059993 


-7.640602 


-10.206973 


-12.766389 


-15.321712 


-17.8 


0.30 


-5.042625 


-7.641730 


-10.221625 


-12.791271 


-15.354531 


-17.9 
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^ neighbor exchange coupling plus the uniaxial single-ion-type anisotropy 
^ energy is assumed for the host-system Hamiltonian. The energy differ- 
^ ; ences between the ground state and the low-lying excited states, which 
^■appear in the Haldane gap, are estimated in the thermodynamic limit. 
^ ! The results are used to analyze the electron-spin-resonance experimental 
data on Ni(C2H8N2)2N02(C104), abbreviated NENP, containing a small 
amount of spin-1/2 Cu 2+ impurities. It is found that in this system, the 
impurity-host (Cu 2+ -Ni 2+ ) coupling is ferromagnetic and its magnitude is 
about 5% of the magnitude of the host-host (Ni 2+ -Ni 2+ ) coupling. 
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§1. Introduction 

In the past years a good deal of work has been devoted to the study 
of the ground and low- lying excited states of the spin-1 antiferromagnetic 
Heisenberg chain. This is mainly motivated by Haldane's prediction 1 ) that 
the ground state of the integer-spin case, in contrast to that of the half- 
integer-spin case, is a massive state characterized by a finite energy gap 
in the excitation spectrum and by an exponential decay of the two-spin 
correlation functions. 2 ) In particular, there has recently been considerable 
interest in the spin- 1/2 degrees of freedom at the edges of the chain with 
open boundary conditions. The spin- 1/2 degrees of freedom are well de- 
scribed by the valence-bond-solid (VBS) picture, 3 ) and lead to the fourfold 
degeneracy of the ground state in the limit of infinite lattice since the ef- 
fective coupling between them vanishes in this limit. This degeneracy has 
been originally found in the so-called AKLT model 3 ) with open bound- 
ary conditions, which is an exactly solvable model including biquadratic as 
well as bilinear exchange interactions. Later, Kennedy 4 ) has shown that 
the ground state of the usual open Heisenberg chain with only bilinear in- 
teractions has the same property. The present authors and Harada 5 ) have 
investigated the low-lying states of the Heisenberg chain with an impurity 
bond, which plays a mediating role between the open and periodic chains, 
in terms of domain- wall excitations. They have shown that a coupling 
between the spin-1/2 degrees of freedom at the edges via the antiferromag- 
netic impurity bond brings about the massive triplet mode in the Haldane 
gap and that this massive mode merges the bottom of the energy contin- 
uum in the limit of the periodic chain. Recently, an appearance of the 
massive mode in the Haldane gap has also been discussed by S^rensen and 
Affleck. 6 ) 
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Hagiwara, Katsumata, Affleck, Halperin, and Renard 7-9 ) have per- 
formed the electron-spin-resonance (ESR) experiment on Ni(C2HgN2)2- 
N02(C104), abbreviated NENP, containing a small amount of spin-1/2 
Cu 2+ impurities, and have given the first experimental evidence for the 
existence of the spin-1/2 degrees of freedom at the sites neighboring the 
impurity spin. They have successfully analyzed their experimental results 
by using the phenomenological Hamiltonian based on the VBS picture, 3 ) 
being composed of three spin-1/2 spins, one of which is the Cu 2+ spin 
and the other two are the spins corresponding to the spin-1/2 degrees of 
freedom. However, they have not discussed the origin of the anisotropy of 
effective exchange interactions between the former and latter spins, which 
they have assumed in their analysis. The spin-1/2 degrees of freedom at 
the edges of the open spin-1 antiferromagnetic chain have also been ob- 
served in NENP doped with nonmagnetic ions such as Zn 2+ , Cd 2+ , and 

Hg 2+ } 10,11) and alg0 

even in 'pure' NENP in which one may find finite mag- 
netic chains broken from infinite chains by, for example, crystal defects. 12 - > 
It has recently been shown 13 -* that the low-temperature behavior of the heat 
capacities and susceptibilities in the NENP:Cu 2+ system is well described 
by Hagiwara et al.'s phenomenological Hamiltonian. 

On the other hand, using a quantum Monte Carlo method, Miyashita 
and Yamamoto 14 -' have studied the open spin-1 antiferromagnetic chain. 
They have calculated the magnetic moment at each site for the lowest 
energy state within the M oc = 1 subspace, M oc being the z-component of 
the total spin for the open chain. According to their results, the one half of 
the total magnetization M oc = 1 is localized in each edge of the chain, and 
the magnitude of the magnetic moment, which takes a maximum value 
at each edge, decays exponentially towards the center of the chain with 

3 



the decay constant of about 6 lattice spacings. Similar results have also 
been obtained by White 15 ) by the density-matrix renormalization-group 
method. Very recently, Yamamoto and Miyashita 16 ) have extended their 
calculation 14 ^ to the case of the presence of the uniaxial single-ion-type 
anisotropy to show that the spin-1/2 degrees of freedom appear also in this 
case, as far as the magnitude of the anisotropy is small enough for the 
ground state to be the Haldane state. 1 ) 

In recent papers 17 ' 18 ) the present authors have analytically investigated 
the equivalence between the above-mentioned Hagiwara et al.'s phenomeno- 
logical Hamiltonian and more realistic Hamiltonians for the NENP:Cu 2+ 
system. Then, they have clearly shown that the single-ion-type anisotropy 
terms in the latter Hamiltonians play an essential role in explaining the 
anisotropy of the effective exchange interactions in the former Hamiltonian. 
This result also confirms the concept of the spin-1/2 degrees of freedom at 
the edges of the open spin-1 chain. 19 ) 

In the present paper, we aim at discussing quantitatively the properties 
of the low- lying excited states as well as the ground state of the spin-1 
antiferromagnetic Heisenberg chain with a spin-1/2 impurity. We express 
the Hamiltonian describing this system as 

H = H oc + H', (la) 

( N-l N } 

n oc = j\ Y^St- s £+1 + dY(si) 2 \ (J>o), (ib) 
I t=i t=i J 

TL' = J' (so • Si + s • Sjv) , (lc) 

where so is the spin-1/2 operator of the impurity (Cu 2+ ) spin; Si is the 
spin-1 operator of the host (Ni 2+ ) spin; TV is the number of host spins; J and 
J 7 are, respectively, the host-host and impurity-host exchange constants; d 
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is the uniaxial single-ion-type anisotropy constant in the host system. Thus, 
Hoc and H' represent, respectively, the Hamiltonian for the host system 
with open boundary conditions and that for the impurity-host exchange 
coupling. It is noted that the Hamiltonian H represents a special (isotropic 
exchange coupling) case of the Hamiltonian which we have discussed in 
ref. 17 as the realistic one for the NENP:Cu 2+ system. 

For the above purpose we employ a method of numerical diagonaliza- 
tion by the Lanczos technique. We calculate, for various values of d and 
J' I J, the energies of the ground and low-lying excited states for finite- 
iV(=5, 7, • • •, 17) systems. Then, we extrapolate these finite- AT results to 
the limit of N — ► oo. Using the results of the extrapolation, we analyze 
the ESR experimental data on the NENP:Cu 2+ system 8 ' 9 ) to determine 
the sign and magnitude of the impurity-host exchange constant J' in this 
system. 

In the next section (§2) the results of the numerical calculation are 
presented and discussed. The analysis of the experimental data is made in 
§3. Section 4 is devoted to a summary. In the Appendix, the expressions 
for the energies of the low-lying states of the Hamiltonian 7i which we have 
obtained in ref. 17 are summarized. We use these in the discussions in 82. 
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§2. Results of the Numerical Calculation 

In order to calculate the energies of the ground and low-lying excited 
states of the system, we first diagonalize numerically, by means of the 
Lanczos technique, the Hamiltonian H for finite N(= 5, 7, 17) sys- 
tems within the subspace determined by the value M = SQ-\-^2d =1 Sf • I n the 
diagonalization we employ the computer program package KOBEPACK/I 
Version 1.0 coded by one of the present authors (M. K.). (The algorithm 
used in this program package is discussed in ref. 20.) Then, we extrapolate 
the finite- TV results to the limit of N — > oo. As shown later, the finite- N 
results up to N = 11 are sufficient to estimate satisfactorily well the limiting 
(N — > oo) values of the energies at least for the parameter region which we 
consider. 

In the present numerical calculation, we choose the following five val- 
ues for d: d = 0.0, 0.1, 0.2, 0.3, and 0.4. For these cfs, the ground state of 
the Hamiltonian H oc is the Haldane state. 21 ) We consider only the ground 
state and the low-lying excited states which appear in the Haldane gap. As 
can be seen from the analytical result obtained in ref. 17, these states are 
those with the lowest three energy eigenvalues within the M — 1/2 sub- 
space and that with the lowest energy eigenvalue within the M = 3/2 sub- 
space. In the following the four eigenvalues for a given value of N are repre- 
sented by £(m,p)(A0, where the different states are denoted by (M, P) with 
P{— +, — ) which stands for the parity with respect to reflection around 
the impurity spin; we have two different (1/2, +) states and the prime is 
added to the state with the larger eigenvalue to distinguish between these 
two states. 

According to our calculation for AT = 5, 7, • • •, 17, the (1/2,+) state is 
always the ground state of the finite chain. We define the energy difference 
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between this state and the other three states as 



A ( m,p) (AO = e (M ,P) (AO - e ( i/2 >+ ) (N) , (2) 

where (M, P) = (1/2,-), (1/2,+)', and (3/2,+). Estimating the limiting 
value A( M p)(oo) of A(m,p)(N), we make a least-squares fit of A( M) p)(5), 
A( M) p)(7), • • •, A( M) p)(17) to the following formula: 

A ( m,p) (N) = A (M;P) (oo) + A exp(-iV/0 W , (3) 

where A, £, and p are constants which take different values depending on 
the three (M, P) states. This procedure leads to the fruitful limiting results 
when 0.02<|J'/J| <0.20 for d = 0.0 and when 0.04 < \J'/J\ <0.20 for 
d>0.0. 

Using A(m,p)(oo), we also define, for convenience, the following quanti- 
ties: 

X / A(i /2 ,_)(oo) 

A(1/2 '- )_ lA (1/2) _)(oo)-A (1/2)+y (oc) 

A , x _ / A (i/2,+)'(oo) 

A _ / A (3/2,+)(00) 

A(3/2, +) - |A (3/2>+ )(00)-A (1/2 , +y (00) 

It should be noted that A^/2,-), A(i/2,+)' and A( 3 / 2 ,+) are the quanti- 
ties which correspond, respectively, to A s (±l/2, 1/2), A t (±l/2, 3/2), and 
A t (±3/2, 3/2) discussed in the Appendix. 

2.1 The case of d = 0.0 

Let us first discuss the case of d = 0.0. Our calculation shows that in 
this case, e ( i /2)+ )/(iV) = e (3/2 , + )(A0 (i.e., A {1/2 ,+y(N) = A (3/2;+) (7V)) holds 
when J'>0 and e( 3 / 2 ,+)(A0 = e(i/2,+)(A0 (i- e -> ^(3/2,+)(N) = 0) holds when 
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for J' > , 
for J' < , 

for J' > , 
for J' < . 

for J' > , 
for J' < , 



(4) 
(5) 
(6) 



J' <0. Thus, the ground state in the J'<0 case is a quartet, in contrast to 
the fact that the ground state in the J'>0 case is a doublet. 

Tables 1(a), 1(b), 1(c) and 1(d) tabulate the finite- N values of 6(1/2, +) (AO, 
e {i/2,-){N) 1 e(i/2,+)' (iV), and e^/ 2j +)(N), calculated for various values of 
J' I J. The values of these energies for N = 17 are plotted as functions of 
J' I J in Fig. 1. We plot the finite-TV values of A (1/2 ,-)(N) and A (1/2 ,+y(N) 
in Figs. 2(a) and 2(b), where the limiting values of these quantities estimated 
by the procedure discussed above are also shown. 

The estimated values of A^/2,-) and &-(\/2,+y (= ^(3/2,+)) are listed 
in Table II, and are plotted versus J' / J in Fig. 3. In order to obtain an 
analytical expression for A( M ,p) in terms of J' /J, we have made a least- 
squares fit of these estimated values to a cubic function 

/j — a (M,P) (J 1 / J) + 0(M,P) (J'IJ? + C( M ,P)(J'IJ? (7) 

with the numerical coefficients cl(m,p), &(m,p) 5 an d c/m,p)- The results, to- 
gether with the value Ahg = 0.4105 J 22 ' 23 ) of the Haldane gap, which is 
defined as the energy difference between the bottom of the energy contin- 
uum and the ground state and is not affected by the presence of an im- 
purity spin, 5 ) are also shown in Fig. 3. This figure should be compared 
with Fig. 2(b) in ref. 17. The coefficients (L(m,p) and 6(m,p) thus determined 
numerically are presented in Table III. As is shown in the Appendix, our 
analytical study 17 ) has led to the results, a(i/2,_) =4/3 and a(\/2,+y = 2, 
which are about 5/4 times as large as the numerical values. As far as the 
ratio a(i/2,-)/«(i/2,+) / is concerned, however, we have a good agreement 
between the analytical and numerical results. Very recently, S^rensen and 
Affleck 6 ^ have obtained the results that a(i/2,-) =ot and ^(1/2,+)' =3o;/2 with 
a~ 1.0640. These are in excellent agreement with our numerical results. 
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2.2 The case of d> 0.0 

As mentioned above, we have carried out in this case the numerical 
calculation for d = 0.1, 0.2, 0.3, and 0.4. In Tables IV(a), IV(b), IV(c) and 
IV(a) we list the finite- AT values of en / 2j+ \(N), en / 2j _\(N), en/2,+y(N), 
and 6(3/2,+) (N) for d = 0.2, calculated for representative values of J' /J. 
The values of these energies for N = 17 are plotted as functions of J' / J in 
Fig. 4. The fmite-AT values of A (1/2 _)(#), A (1/2y+y (N), and A (3/2>+) (JV) 
for d = 0.2 as well as the limiting values of these quantities estimated by the 
above procedure are plotted in Figs. 5(a), 5(b) and 5(c). 

The estimated values of An /2,-)> ^(i/2,+)'i and A( 3 / 2 ,+) for the above- 
mentioned values of d are tabulated in Table II, and those for d= 0.2 are 
plotted versus J' / J in Fig. 6. In the present case, we have also obtained 
the analytical expression for A( M p) by the same method as that used in 
the case of d = 0.0, and the results for the coefficients cl(m,p) and 6(m,p) 
are listed in Table III. As we will discuss in §3, the Haldane gap Ahg for 
d=0.2 is given by Ahg = 0.30 J. This value of Ahg as well as the analytical 
expression for A( M; p) obtained for <i = 0.2 is shown in Fig. 6. This figure 
should be compared with Fig. 2(a) in ref. 17. 
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§3. Analysis of the Experimental Data 

We now turn to the analysis of the ESR experimental data on the 
NENP:Cu 2+ system. 8 ' 9 ) First we discuss the value of d for NENP. Go- 
linelli, Jolicceur, and Lacaze 24 ) have numerically analyzed the splitting of 
the triplet state, which lies at the bottom of the energy continuum in the 
case of the periodic chain, into singlet and doublet states due to the uniaxial 
single-ion-type anisotropy term d^2^ =1 (S%) 2 . According to their result, the 
excitation energy A s of the singlet state measured from the singlet ground- 
state energy and the corresponding excitation energy Ad of the doublet 
state are given by 

A s = (0.41 + 1.41d)j, (8) 
A d = (0.41 -0.57d)J, (9) 

for 0.0 < d<0.3. Golinelli et al. 24 ^ have compared the inelastic-neutron- 
scattering (INS) experimental results 25 ) A s = 2.5 meV and Ad = 1.15 meV 26 ) 
on NENP with eqs. (8) and (9) to obtain 

J = 3.75meV = 43.5K and d = 0.18. (10) 

On the other hand, if we use the results of magnetization measurements, 27 ^ 
A S = 27.5K and A d = 11.4K, 28 ) we obtain 

J = 39.1K and d = 0.21 . (11) 

Thus, we adopt here d=0.2 as the value of d for NENP. 

Hagiwara et al.'s ESR experimental results 8 ' 9 ^ show that the ground 
state and the first, second, and third excited states in the NENP:Cu 2+ 
system are the M = ±1/2, ±3/2, ±1/2, and ±1/2 states, respectively. 
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From this together with our analytical results 17 ^ (see Fig. 2 in ref. 17 and 
eqs. (A-l), (A-2), (A-3) and (A-4) in the Appendix), we can conclude that J' 
is ferromagnetic (J 7 <0). The energy separations between the three excited 
states and the ground state are A^ xp = 5.7GHz, A 2 xp = 25.8GHz, and Ag Xp = 
71.7GHz. These correspond, respectively, to A( 3 / 2) +)(oo), A( 1 / 2) _)(oo), and 
^(1/2,+)' (oo). From Table III and eqs. (4-6) we have, for -1<J'/J<0, 

A( 3 /2,+) (oo)/ J = -0.192 ( J' I J) - 0.503 (J 7 / J) 2 , (12) 
A(i /2 _)(oo)/J = -0.592 (J' I J) + 0.125 (J'/J) 2 , (13) 
A (1/2 ,+)/(oo)/J= -1.618 (J'/J) - 0.293 (J'/J) 2 . (14) 

On the other hand, the Haldane gap Ahg for finite d is given by Ad, since 
it should be defined as the energy difference between the bottom of the 
energy continuum and the ground state. 5 ) Thus, Ahg = 0.30 J fc> r d=0.2, 
while its experimental value we use here is A^q = 1.15 meV = 278 GHz. 25 ) 
We determine the value of J'/J from 

A (1/2 , +) ,(oo):A H G = Ar P :ASS, (15) 

which gives 

J'/J =-0.048. (16) 

Then, we have 

A( 3 /2,+)(oo) : A (1/2 _)(oo) : A (1/2;+) , (oo) : A HG 

= 0.03 : 0.10 : 0.26 : 1.00, (17) 

and this ratio is in very good agreement with 

A exp . A exp . A exp . = q q 2 . Q.Q9 I 0.26 : 1.00 . (18) 
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§4. Summary 

We have investigated the ground and low-lying excited states of the 
spin-1 antiferromagnetic Heisenberg chain with a spin-1/2 impurity, the 
Hamiltonian of which is given by eqs. (la-c). Employing the computer pro- 
gram package KOBEPACK/I Version l.O, 20 ^ we have calculated, for various 
values of d and J' /J, the energies of these states for finite- N(= 5, 7, • • •, 
17) systems. We have extrapolated to the limit of N — > oo the finite-TV 
results for the energy difference between the ground state and the low-lying 
excited states, assuming the formula given by eq. (3), to estimate the values 
in this limit. Using these limiting results, we have analyzed Hagiwara et 
a/.'s ESR experimental data on the NENP:Cu 2+ system. 8 ' 9 ) We have found 
that in this system, the magnitude of the impurity- host (Cu 2+ -Ni 2+ ) ex- 
change constant is 0.048 times as large as that of the host-host (Ni 2+ -Ni 2+ ) 
exchange constant, the former being ferromagnetic. 
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Appendix 

In ref. 17 we have calculated the energies of four low-lying states of H, 
confining ourselves to the Haldane region where the ground state of 7i oc is 
the Haldane state, and treating TV as a small perturbation (| J'\j J<C l) . 29 - ) 
In the calculation, we have employed, as the unperturbed wave functions, 
the fourfold-degenerate ground state wave functions of Ti oc obtained in our 
previous paper, 5 ) which are expressed in the matrix-product form. 30 ' 31 ) The 
calculated energies, measured from the ground state energy Eo? (see eq. (7) 
in ref. 17) of Ti oc , are given by 

/ 3 3\ T/ 4-d 1X 




1 + 16§±^-1>, (A- 2) 



4^) =o (A - 4) 

in the Haldane region 4 > d > — 2. Here, we have denoted the energy of 
the state belonging to the subspace determined by M as e r (M, 5) (r = s, 
t), where S represents the magnitude of the total spin of the corresponding 
state in the isotropic (d = 0) case, and the subscripts s and t show that 
the energies are associated with the singlet and triplet states for Ti oc , re- 
spectively If we choose the origin of the energies ^0? + £t(±l/2, 1/2) and 
define A r (M, S) (r = s, t) as 

A r (M, S) = s v (M, S) - e t (±\, \y (A • 5) 

then An/2,-), ^(1/2,+y an d ^(3/2,+) correspond, respectively, to A s (±l/2, 
1/2), A t (±l/2, 3/2), and A t (±3/2, 3/2). 
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Figure Captions 

Fig. 1. Plots versus J' / J of 6(m,p)(17) for d = 0.0. The open circles 
show the numerical results for 6(1/2, +) (N), the crosses show those for 
e (i/2,-){N) 1 the open triangles show those for e(i/2,+y(N), and the 
solid circles show those for e( 3 / 2 ,+) (N). The solid lines are the guide 
to the eye. 

Fig. 2. Plots versus J' /J of A (M)P) (iV) for d = 0.0: (a) A (1/2 _)(iV) and 
(b) An /2,+y (N). The open circles show the numerical results for finite- 
N systems; iV = 5, 7, • • •, 17 from the topmost circle to the lowest circle 
for each value of J' / J. The solid lines are the guide to the eye. The solid 
squares show the limiting (TV^oc) results estimated by the procedure 
discussed in the test. 

Fig. 3. Plots versus J' / J of A(m,p)(N) for d = 0.0. The crosses show 
the estimated values of A(!/ 2 ,-) and the open triangles show those of 
An /2,+y (= ^(3/2,+))- The solid lines are least-squares fits to eq. (7). 
The dashed line represents the bottom of the energy continuum. 

Fig. 4. Plots versus J' / J of €(m,p)(17) f° r d = 0.2. The open circles 
show the numerical results for en / 2) +)(iV), the crosses show those for 
e (i/2,-){N) 1 the open triangles show those for e(i/2,+y(N), and the 
solid circles show those for £(3/2,+) (N). The solid lines are the guide 
to the eye. 

Fig. 5. Plots versus J' /J of A (M , P) {N) for d = 0.2: (a) A (1/2 _)(iV), 
(b) A( 1 / 2)+ )/(AT), and (c) A( 3 / 2)+ )(iV). The open circles show the nu- 
merical results for finite- N systems; N — 5, 7, • • •, 17 from the topmost 
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circle to the lowest circle for each value of J' / J. The solid lines are the 
guide to the eye. The solid squares show the limiting (TV^oo) results 
estimated by the procedure discussed in the test. 

Fig. 6. Plots versus J' / J of A( M)P )(iV) for d = 0.2. The crosses show the 
estimated values of An/2,-), the open triangles show those of An/ 2) +y , 
and the solid circles show those of A( 3 / 2) +). The solid lines are least- 
squares fits to eq. (7). The dashed line represents the bottom of the 
energy continuum. 
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